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We study the ground state properties of rectangular quantum dots by using the spin-density- 
functional theory and quantum Monte Carlo methods. The dot geometry is determined by an 
infinite hard-wall potential to enable comparison to manufactured, rectangular-shaped quantum 
dots. We show that the electronic structure is very sensitive to the deformation, and at realistic 
sizes the non-interacting picture determines the general behavior. However, close to the degenerate 
points where Hund's rule applies, we find spin-density-wave-like solutions bracketing the partially 
polarized states. In the quasi-one-dimensional limit we find permanent charge-density waves, and 
at a sufficiently large deformation or low density, there are strongly localized stable states with a 
broken spin-symmetry. 

PACS numbers: 73.21. La, 71.10.-w 



I. INTRODUCTION 

During the rapid development of nanotcchnology, the 
discoveries in the physics of small electronic structures 
have concurrently opened new channels in this extremely 
active field of both theoretical and experimental research. 
Quantum dots, which fundamentally are confined elec- 
tron bunches, represent basic components of nanoelec- 
tronics. They have been shown to possess many atom- 
like properties, such as the specific shell structure, deter- 
mined by the properties of the external confinement. 1 

In lithographically fabricated quantum dots the elec- 
trons are strictly confined on the interface of the semicon- 
ductor heterostructure, which makes the dot essentially 
two-dimensional. The lateral confinement, created by 
adding a voltage to the top gate electrodes, is most com- 
monly approximated by the harmonic oscillator poten- 
tial. In the isotropic case, this modeling has been shown 
to lead to a similar addition energy spectrum as mea- 
sured in the experiments, and adjustments in the model 
potential have made the agreement even more precise (see 
Ref. 2 for a review). 

Deviations from parabolic confinement have most com- 
monly been studied in connection with the far-infrared 
response (FIR)i^^iSii This is due to the generalized 
Kohn's theorem(2*S stating that FIR couples only to the 
center-of-mass (CM) motion which in the case of a perfect 
parabolic potential can be separated from the relative 
motion (RM). Since the CM motion has the same energy 
eigenvalues and dipole resonance frequencies as a single 
electron, no information on internal degrees of freedom 
can be obtained. Ugajin 6 studied FIR absorption for a 
two-electron square-well quantum dot by using exact di- 
agonalization, and recent density-functional calculations 
of corner and side modes for triangular and square dots 
have been done by Valfn-Rodriquez et ali 

The ground state electronic structure of square-shaped 
quantum dots was first calculated by Bryant^ He used 
configuration-interaction (CI) methods to examine the 



role of interactions for two electrons. Creffield et alH* 
studied polygonal two-electron quantum dots with nu- 
merically exact diagonalization, concentrating on the 
Wigner crystallization^ i.e., the localization of the elec- 
trons due to the dominant Coulomb interaction in the 
low-density limit. In our previous work^ we found an 
agreement with their results for polygonal dots by using 
the spin-density- functional theory (SDFT). We extended 
the examination to larger electron numbers, including 
broken spin-symmetry configurations. Those states cor- 
respond to spin-density waves (SDW) found in the weak- 
confinement limit of parabolic quantum dots and rep- 
resent energetically stable and accurate solutionsi 2 *^*^ 
Akbar and Le©i& also used the SDFT to calculate the 
addition energy spectrum for square quantum dots with 
different sizes. 

Until now, the study of square-shaped quantum dots 
with a hard-wall confinement has not been generalized 
into arbitrary rectangular shapes. However, experiments 
have been done on rectangular mesas of vertical dots 
by Austing et ai.pl who applied electron-beam lithog- 
raphy with etching techniques on a double barrier het- 
erostructure (DBH)>i^ They measured the addition spec- 
trum with different deformation parameters as well as the 
magnetic field dependence on the Coulomb oscillations. 
In the same extensive study, they performed SDFT cal- 
culations to simulate the external confinement with an 
elliptic potential. That approximation was shown to be 
tentative, though insufficient for a general description of 
rectangular quantum dots. Lee et ali^ also studied ellip- 
tical dots with the SDFT, including additional harmonic 
confinement in the z direction, and obtained similar ad- 
dition energy spectra. 

In the present paper, our secondary aim is to test the 
ability of a hard-wall external confinement to approx- 
imate real rectangular quantum dots, measured in the 
above-mentioned study. Our main purpose is, however, 
to clarify the electronic behavior in a rectangular box, 
beginning from a basic textbook-example of quantum 
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mechanics and leading to the discussion of the role of 
interactions and symmetry-broken solutions in different 
regimes. 

The outline of the paper is as follows. In Sec. ITT1 we 
present the model Hamiltonian and the analytical shell 
structure of a two-dimensional rectangular box. The 
computational methods, a real-space SDFT technique 
and the variational Monte Carlo (VMC) method are in- 
troduced in Sec. IIIII From the results in Sec. IIVI we 
first give the chemical potentials and the addition energy 
spectra of rectangular quantum dots. Then we continue 
toward a deeper insight into the electronic structure, in- 
cluding the spin-behavior in the dot and the quasi-one- 
dimensional limit. The paper is finished with a summary 
in Sec. 
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II. MODEL AND THE SHELL STRUCTURE 

We define our quantum dot to be two-dimensional, i.e., 
strictly confined in the z direction. We use the effective 
mass approximation (EMA) to describe electrons mov- 
ing in the plane, surrounded by background material of 
GaAs with the effective electron mass m* = 0.067m e and 
dielectric constant e = 12.4. Energies are thus given in 
Ha* w 11.8572 meV and lengths in k 9.79 nm. 

The model Hamiltonian of an A-electron system in an 
external potential can be written as 
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The external confinement in the xy plane is described by 
an infinite hard-wall potential, 



V ext (x,y) 
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oo, elsewhere. 



(2) 



Therefore, the area of the dot is PL 2 , where the defor- 
mation parameter (3 defines the ratio between the side 
lengths of the rectangle. 

Let us now omit the mutual interactions of the elec- 
trons, and consider the single-electron states in a two- 
dimensional rectangular box. We need two quantum 
numbers, n x and n y , to label all the needed eigenfunc- 
tions of two Cartesian coordinates. Inside the box we can 
write an explicit formula for these functions as 
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Inserting the eigenfunctions to the stationary 
Schrodinger equation and setting the area of the 
rectangle PL 2 = it 2 , give now the energy eigenvalues in 
a simple form, 



En x ,n y = U^+Pnl 



(4) 



FIG. 1: Lowest single-electron eigenenergies for rectangular 
quantum dots as a function of the deformation. 
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FIG. 2: Sum of the lowest energy eigenvalues, V" e„ „ 



for N 



10, 12, and 14 non-interacting electrons as a 



function of the deformation. 



Fig- HI shows these eigenvalues as a function of p. The 
degeneracies in the case of j3 = 1 introduce the magic 
electron numbers for a square, N = 2,6, 8, 12, 16, 20, . . ., 
corresponding to closed shells. When the dot is squeezed, 
the degeneracies are lifted, resembling the behavior of 
the single-electron states in an anisotropic harmonic os- 
cillator potential. In the rectangular case, however, one 
cannot find such regularly located junctions of the eigen- 
states as in elliptic dots. This is a direct consequence 
of the more constricted symmetry of rectangular than 
harmonic quantum dots. This produces remarkable dif- 
ferences in the electron structures as will be shown be- 
low. Fig. |21 gives the sums of the lowest eigenvalues for 
N = 4, 6, 14. We can find formation of stable configu- 
rations with certain (N, P) combinations as local minima 
in the total energy curve. Correspondingly, the cusps 
indicate degeneracies of the states. 

Accumulation of states (n x ,l) at high deformation, 
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which can be seen in Fig. ^ is similar to the forma- 
tion of Landau bands in the Fock-Darwin energy spectra 
for the harmonic oscillator potential at high magnetic 
fields^ As the deformation is made stronger, the system 
becomes gradually quasi-one-dimensional and the occu- 
pation of the electrons is determined by the quantization 
in the longer direction. 



III. COMPUTATIONAL METHODS 
A. Spin-density-functional theory 

We employ the usual self-consistent formulation of 
the density-functional theory, introduced by Kohn and 
Sham^i The single-electron wave functions are solved 
within the EMA from 
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where the effective potential is a sum of the external, 
Hartree, and exchange-correlation potentials, 



KffM = Kxt(r) + Mr) + V£(r). 



(6) 



To calculate V° c (r), we use the local spin-density ap- 
proximation (LSDA), 

6E LSDA f 
V x l(r) * -j^y = J drn(r)e xc (n(r),C(r)), (7) 

where e xc is the exchange-correlation energy per electron 
in a uniform electron gas of density n = n-f + n± and 
spin-polarization £ = (n-f — ni)/n. We employ a recent 
analytic paramctrization for e xc , formulated in connec- 
tion with diffusion Monte Carlo calculations (DMC) by 
Attaccalite et al*& It is written as 

e xc (r s ,C) = e x (r s ,C) + (e~^ - l)ef(r s ,C) 

+ a {r s )+a 1 {r s )( 2 +a 2 (r s )(\ (8) 

where r s = l/y/irn is the density parameter for the 
2DEG, a's are density dependent functions of the gen- 
eralized Perdew-Wang formjSi (3 = 1.3386, and e x is the 
exchange energy given as 

e x (r s , C) = -2V2[(1 + C) 3/2 + (1 - C? /2 ]/^r s . (9) 

In Eq. (JHJ , e x 6 ' is the Taylor expansion of e x beyond the 
fourth order in £ at £ = 0. 

The above parametrization fits to the DMC simula- 
tions over the whole range of spin-polarization (0 < 
C < 1). This is an essential extension to the often- 
used parametrization of Tanatar and Ceperleyj 2 * which 
is based on DMC calculations for systems with £ = 
and 1. Gori-Giorgi et alm^ have shown that the improve- 
ment gained with the new parametrization is directly pro- 
portional to the electron density and the polarization of 



the system. In our recent article^ we compare different 
LSD functionals in small 2D quantum dots. We show 
that in comparison with the VMC calculations, the new 
parametrization by Attaccalite et al. gives more accu- 
rate results for the exchange-correlation than the forms 
of Tanatar and Ceperley. 

We perform the numerical calculations in real space 
with two-dimensional point grids without implicit sym- 
metry restrictions. Through this approach, we can shape 
the external potential almost arbitrarily in the comput- 
ing region. The number of grid points is 128 x 128, which 
gives an accuracy of better than ~ 1% in the total en- 
ergy. To accelerate the numerical process, we apply the 
Rayleigh Quotient Multigrid (RQMG) method 2 ^ for the 
discretized single-electron Schrodinger equation (jSJ). A 
detailed description of this method, generalized to an ar- 
bitrary number of lowest eigenenergy states, can be found 
inRef.lM 



B. Variational quantum Monte Carlo method 

The variational quantum Monte Carlo (VMC) 29 
method starts from constructing a trial many-body wave 
function '5 with desired properties and with free varia- 
tional parameters a;. The parameters are then optimized 
to converge toward the exact wave function 'J'o- Using 
the optimized wave function, the expectation value of an 
observable A can be evaluated as the average of the cor- 
responding local quantity i f?~ 1 A i Sf. For example, energy 
is found from the Hamiltonian operator H as: 



E* = lim ^y^M = (*|H|*) 



(10) 



where the V-particlc- coordinate configurations Ri are 
distributed as |\l/| 2 and generated using the Metropolis 
algorithm. 

The variational principle guarantees that the total en- 
ergy given by the VMC method, using any trial wave 
function with proper particle symmetry, is always an up- 
per bound for the true total energy of the quantum state 
in question. The variance of the local energy ^~ 1 H'^> di- 
minishes as the trial wave function approaches an eigen- 
state of the Hamiltonian, and as a result it can be used 
not only as a measure of the statistical error in E^, , but 
also as a measure of the difference between calculated 
and true energies E^ — E^ . 

The variational parameters in the trial wave function 
are optimized by minimizing the total energy. The mini- 
mization process itself was done using the stochastic gra- 
dient method^ The method has proven to be fast and 
reliable. 

The variational wave functions used in this work are 
of the form 
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where the two first factors are Slater determinants for 
the two spin types, and J is a Jastrow two-body cor- 
relation factor. We neglect the three-body and higher 
correlations. This has shown to be very accurate in our 
previous VMC studies (See, e.g., Refs. l3lll33l3^ . For the 
Jastrow factor we use 
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where a is fixed by the cusp condition to be 3 for a pair of 
equal spins and 1 for opposite ones and b is a parameter, 
different for both spin-pair possibilities. C is the scaled 
Coulomb strength. The single-particle states in the de- 
terminants are taken to be those for the non-interacting 
problem given in Eq. ©. 



IV. RESULTS 
A. Addition energy spectra 

We calculate the total energies of rectangular dots 
with different deformation parameters up to 16 elec- 
trons. We keep the dot area constant, A = it 2 , through 
our calculations. The density parameter, defined as 
r s = y/A/(Nn) (Ref.E3), thus gets values between 0.44 
and 1.8. The electron density in our quantum dots is 
therefore higher on the average than that of Austing et 
alH with r s — 1.5. Nevertheless, we find that the differ- 
ence has no noticeable effect on our results. In Fig. [3] we 
show the SDFT and VMC results for the chemical poten- 
tials, fi(N) = E(N) - E(N - 1), (E(0) set to zero), for 
various values of ft. The agreement between the results 
is good and independent of N and ft. As a consequence 
of the two- fold degeneracy in the eigenstates, the pair- 
ing of chemical potential values dominates the picture. 
However, a closer look reveals deviations from this ten- 
dency. Due to Hund's rule, near the degenerate points in 
the single-electron spectrum, the spins of the two highest 
energy electrons are parallel and they occupy different 
states. So, there are regimes in which fi(N + 1) and 
fi(N — 1) behave in the same way, for example N = 8 
as ft ~ 1.3 — 1.5, corresponding to the degeneracy of the 
states (n x , n y ) = (2,2) and (1,3). Similar effects in chem- 
ical potentials have been observed in measurements of 
vertical quantum dots in magnetic fields^ and in calcu- 
lations of elliptically deformed dots^ Due to the rather 
coarse spacing of our ft- values in Fig- El all the deviations 
are not observable. A more detailed description as well 
as a comparison to elliptic dots follows below. 

In Fig. 0| we show the addition energies, fi(N + 1) — 
fi(N), tor ft — 1 — 4. The spectra obtained with the 
SDFT and VMC coincide well, especially for N < 10. In 
the case of a square dot (ft — 1), the magic configurations 
can be seen as large peaks in the spectrum. For ft = 1, the 
relatively large addition energy for N = 4 corresponds to 
a half-filled shell according to Hund's rule. The spectrum 




FIG. 3: Calculated chemical potentials for rectangular quan- 
tum dots as a function of the deformation parameter. The 
SDFT and VMC results are given by pluses and crosses con- 
nected with solid and dotted lines, respectively. 



agrees well with the results of Akbar and LeeiS for a 
square dot of a similar size. 

In general, the results for the rectangular quantum dot 
are very sensitive to the deformation. As ft increases, the 
peaks for N = 4, 12 rapidly vanish but reform above ft ~ 
1.2. For N = 8, 14 the addition energy oscillates more 
smoothly, and in dots with N = 6, 10 it varies relatively 
slowly, declining in the former and growing in the latter 
between ft = 1 — 1.5. Above ft ~ 2 the formation of 
an even-odd structure corresponds to the filling of states 
(n x ,l). In that regime, the growing amplitude in the 
peaks reflects the increasing spacing between the single- 
electron eigenstates shown in Fig. ^ 

It is intriguing to compare qualitatively the evolution 
of the spectra in the regime of /? ~ 1.3 — 1.5 to the ex- 
perimental results of Austing et al£L There are two diffi- 
culties in the direct comparison. First, the experimental 
mesa is much larger than the area where the electrons 
are actually confined, causing uncertainty of the value for 
the deformation parameter. Secondly, there are evidently 
irregularities in the experimental dots, leading to unex- 
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FIG. 4: Addition energy spectra for rectangular quantum dots 
with different deformation parameters. The SDFT and VMC 
results are given by pluses and crosses connected with solid 
and dotted lines, respectively. 



pected behavior in the spectrum as speculated by Aust- 
ing et alH In spite of these problems, we can generally 
find similarities in the spectra. Compared with the ellip- 
tic case, there is more tendency of forming peaks for even 
N in both the experiments and our approximation at 
(3 ~ 1.3 — 1.5. This may result from the higher symme- 
try of the elliptic than rectangular dot, discussed in the 
context of the single-electron spectrum in Sec.[n] Of par- 
ticular electron numbers, the behavior of the curves for 
N = 2, 6, and 10 qualitatively agree, and the biggest dif- 
ference is the rapid disappearance of the peak for N = 4 
at f3 = 1.44 — 1.5 in the experiment. Overall, our hard- 
wall approximation seems to be a slightly better approx- 
imation for rectangular-shaped quantum dots than the 
elliptic description in Ref. [13 However, more accurate 
comparison than presented here would certainly require 
more measurements and over a wider range for f3. An 
ideal experimental setup would also contain a way to 
tune (3 for a single dot, reducing the variation induced 
by using different dots for different (3. 



B. Spin-development and the role of interactions 

Next we consider more carefully the effect of electron- 
electron interactions on the electronic structure. In Fig-El 
we compare the non-interacting single-electron spectrum 
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FIG. 5: Non-interacting eigenenergies and the ground state 
spins of N interacting electrons for rectangular dots as a func- 
tion of the deformation f3. 



with the evolution of the total spin for even N . Due to 
Hund's rule, we can see partial spin-polarization (5 = 1) 
close to every degenerate point in the single-electron en- 
ergy spectrum. In the case of a square, the S = 1 
ground state is found correctly for half-filled shells with 
N = 4, 10, and 14. The spin state changes rapidly to 
5 = as the dot is squeezed. The range of 5 = 1 
regimes is obviously directly proportional to the slope dif- 
ferences of the crossing eigenstates. The triple-crossing 
for N = 16 at (3 « 1.7 leads to two separate 5 = 1 
regimes around the degenerate point. In most cases, po- 
larization occurs at higher (3- values than the correspond- 
ing crossing of the non-interacting states. Therefore, by 
taking the electron-electron interaction into account, the 
effective deformation of the rectangle is lower than that of 
the bare external potential. This is contrary to the result 
for elliptic dots obtained by Lee et al.^ who concluded 
that the interactions tend to strengthen the bare poten- 
tial by a factor of ~ 1.15 — 1.25. Intuitively, one would 
expect just an opposite behavior: in hard-wall rectangu- 
lar dots the maximum electron density is pushed toward 
the shorter sides, whereas elliptic and harmonic confine- 
ments favor pronounced density at the center. We will 
present this tendency explicitly in Sec. IIV CI 

As we show in Fig. every 5=1 state is bracketed 
by spin-density-wave-like solutions. In these regimes, the 
exchange-energy gained in the polarized state is relatively 
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FIG. 6: Spin-polarization of a 12-electron rectangular quan- 
tum dot in two SDW regimes. 
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The phenomenon has an analogy in molecular systems, 
known as the spontaneous Jahn- Teller effects^ 5 , any non- 
linear molecular system in a degenerate electronic state 
will be unstable and will undergo distortion to form a sys- 
tem of lower symmetry and lower energy. In this particu- 
lar case, however, the commonly used argument that the 
symmetry-broken state would make the electronic struc- 
ture more stable by the enlargement of the Fermi gap is 
not precisely valid, as can be concluded from Fig. It 
is more or less a matter of preserving the S = state 
against the transition to the S — 1 state, representing 
here a more stable configuration. 

The SDW state is a mixture of different S* = states, 
and there has been a lot of debate if this mixed state is 
physically meaningful. 2 In our forthcoming studies which 
include the exact diagonalization results for N = 4, we 
hope to enlighten the validity of the above-represented 
mixed states for rectangular quantum dots. Until now, 
however, symmetry-broken solutions have shown their el- 
igibility in several systems. In parabolic quantum dots, 
for example, the SDW state was found to agree aston- 
ishingly well with VMC results in the weak-confinement 
limit, 15 especially when the latest 2D-LSDA functional 
was usedi^i In our previous study we showed that 
in polygonal quantum dots the breaking of the spin- 
symmetry precedes the complete Wigner molecule for- 
mation at low densities^ Quantum wires, studied in the 
context of SDW solutions by Reimann et al.j& represent 
another interesting example that we discuss in the next 
section. 



FIG. 7: Development of the Kohn-Sham energy levels for a 
rectangular eight-electron dot as a function of /3. Doubly 
degenerate levels are denoted by (2) and the occupied states 
by (o). 



close to the cost paid by occupying the higher energy 
state. By breaking the internal spin-symmetry, the dot 
gains exchange-correlation energy which preserves it at 
the paramagnetic 5* = state instead of following Hund's 
rule with S = 1. A similar behavior was found in the 
study of elliptic dots for certain configurations. 17 In the 
resulting SDW-like solution, the spin-up and spin-down 
densities are symmetrically coupled with each other as 
shown in Fig. [fj] for a 12-electron dot with j3 — 1.14 and 
1.8, corresponding to two symmetry-broken regimes. In 
both cases there are six maxima and six minima in the 
spin-polarization but the shapes of the waves are totally 
different. 

Besides electron densities, it is interesting to consider 
the development of the Kohn-Sham energy levels near 
the degenerate point. In Fig. 0we show the evolution 
for an eight-electron dot with (3 — 1.2,1.3,1.4, and 1.5. 
As can be seen in in Fig. these values correspond to 
states S = 0, SDW, S = 1, and SDW, respectively. In 
the SDW states, the Fermi gap is just large enough to 
prevent the polarization on the highest occupied level. 



C. Quasi-one-dimensional limit 

As the deformation is made larger, electrons in the 
dot become gradually restricted in the lowest energy 
state in the y direction, i.e., only states (n Xl 1) are filled. 
This corresponds to the quasi-one-dimensional limit and 
a quantum-wire-like electronic structured Beyond this 
limit, we find two phases directly observable in the elec- 
tronic density. First, there is a charge-density wave 
(CDW) with N/2 peaks and preserved spin-symmetry. 
As the deformation or the dot size is increased further, 
a spin-density wave appears, consisting of interlocked 
spin-up and spin-down contributions and resulting in a 
Wigner-molecule-like electron density with N peaks. In 
Fig. IS] we show examples of both cases with electron 
density profiles and the corresponding Kohn-Sham eigen- 
functions. In both wires, the area is still 7r 2 , correspond- 
ing to r s = 0.51 (N = 12) and 0.63 (N = 8). As can 
be seen in the figure, the 12-electron wire with f3 = 10 
retains the spin-symmetry and the KS eigenfunctions are 
doubly degenerate. On the contrary, the 8-electron wire 
with [3 = 20 has a broken spin-symmetry. In this case 
the single-electron KS eigenfunctions are mirror images 
of each other, and therefore the KS energy levels are still 
doubly degenerate. Due to the dominating Coulomb in- 
teraction, the lowest KS eigenfunctions correspond to lo- 
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FIG. 8: Electron density profiles (a,c), and the corresponding 
single-electron Kohn-Sham eigenfunctions (b,d) for two quan- 
tum wires with spin-symmetry-preserved and broken ground 
states. The parameters (N,/3,r 3 ) = (12, 10,0.51) in (a,b) and 
(8, 20, 0.63) in (c,d). Doubly degenerate functions are denoted 
by the numbers (2). 
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FIG. 9: Phase separation curves between the spin-symmetry- 
preserved and broken solutions for rectangular quantum dots 
in the quasi-one-dimensional limit. 



calized states near the ends, having 0.16 Ha* lower energy 
than the other occupied levels with a mutual separation 
of ~ 0.06 Ha*. Compared to this, the Fermi gap is par- 
ticularly large, 0.69 Ha*. In this sense, the breaking of 
the spin-symmetry resembles two-dimensional systems in 
the low-density limit i& 

Next we vary both the value of (3 as well as the density 
parameter r s in order to examine the transition point be- 
tween the two phases discussed above. As shown in Fig-El 
for N — 4 and 6, the r s value needed for the deformation 
to form a SDW decreases as (3 is increased. The behavior 
is rather insensitive to N; larger electron numbers also 
qualitatively follow the presented curves with the same 



tendency of moving slightly up in r s with increasing elec- 
tron number, which may arise from our definition for r s . 
We point out that beyond the phase separation, i.e., at 
particularly large j3, r s values, there is also a transition to 
a fully spin-polarized state not shown in the figure. For 
example, for a four-electron wire with /3 = 10, this occurs 
at r s ~ 1. In three-dimensional metal nanowires, SDFT 
calculations have similarly been shown to lead to spon- 
taneous polarization in zero magnetic field at a critical 
radius of the wire^ Recent conductance measurements, 
performed for ultra-low-disorder quantum wires, support 
this phenomenon^ Another remark concerning Fig. is 
the fact that the ability to reach the quasi-lD-limit re- 
quires naturally smaller electron number than there are 
available lowest (n x , 1) states to be occupied. This con- 
dition can be easily estimated from the non-interacting 
single-electron spectrum (Fig. ^ . 

Comparison between our rectangular hard-wall quan- 
tum wires and elliptical wires with harmonic confinement 
studied by Reimann et alm& reveals some noticeable dif- 
ferences. First, our LDA (spin-compensated) solution 
always has a CDW with N/2 pronounced maxima, con- 
trary to the elliptic case with a smooth electron density. 
Secondly, in rectangular wires the total density distribu- 
tion is remarkably concentrated at the ends due to the 
dominating lowest KS eigenstates shown in Fig. [S] The 
opposite distribution in these two geometries is a direct 
consequence of the difference in the confining potential: 
in the elliptic wire, the bowl-like restriction along the 
wire accumulates a pronounced density at the center, 
whereas in the hard-wall wire the Coulomb interaction 
pushes the dominant distribution to the ends. Increasing 
f3 or r s emphasizes this tendency of localization. It is no- 
ticeable that the SDW formation is the origin of both the 
particularly large Fermi gap and the strong localization 
of the lowest eigenfunctions. 



V. SUMMARY 

We have investigated the electronic properties of hard- 
wall rectangular quantum dots. Most calculations have 
been performed with a symmetry-unrestricted SDFT 
scheme in real-space. For the addition energy spectra, 
we have done also VMC calculations and found excellent 
agreement between the two methods. Direct comparison 
with experiments for rectangular mesas of vertical quan- 
tum dots is troublesome but we find tentative common 
features in the addition energy spectra. Close to the de- 
generate points where Hund's rule applies, the states with 
partial polarization are bracketed by unstable SDW-like 
solutions. The effective deformation is generally lower 
than that of a bare potential, but the general picture 
follows the non-interacting single-electron spectrum. Be- 
yond the quasi- ID-limit we find very stable SDW states 
and extremely strong localization near the ends of the 
wire, arising from the shape of the hard- wall confinement. 
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